Predicted giant magnetic moment on non-{nOm} surfaces of d-wave superconductors 

— Can it be observed and how? 



Chia-Ren Hu 1 and Xin-Zhong Yan 1,2 
1 Department of Physics, Texas A&M University, College Station, TX 77843-4242 
2 Institute of Physics, Chinese Academy of Sciences, P. O. Box 603, Beijing 100080, China 

(February 1, 2008) 

It has been noted previously that the sizable areal density of midgap states which must exist on 
any non-{n0m} surface of a d-wave superconductor can lead to a giant magnetic moment. Here 
we show that this effect is observable, and discuss two precise ways to observe it: (i) by directly 
measuring magnetic moment in a system with a large density of internal {110} surfaces, or (ii) by 
performing spin-polarized tunneling on a {110} surface. In both cases a sufficiantly large magnetic 
field should be applied in the [110] direction. 
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One of us (CRH) has noted previously [1] that the 
sizable areal density of "midgap states" (MS's) which 
must exist on any non-{nOm} surface of a d-wave super- 
conductor (DWSC) (with n and m integers or zero) can 
lead to a "giant magnetic moment" (GMM). The MS's 
are topological signatures of unconventional pairing sym- 
metry. They are "nearly-dispersionless" quasi-particle 
states, characterized by momenta along the surface rang- 
ing from — kp to Uf (&f being the Fermi momemtum), 
but all with the same "zero" energy as measured from 
the Fermi energy (in the WKBJ approximation). These 
states can lead to a narrow dcnsity-of-states (DOS) peak 
at the Fermi energy, where the integrated bulk-DOS dips 
to zero. One of the observable consequences of these MS's 
is therefore a "zero bias conductance peak" (ZBCP) in 
single-particle tunneling [1,2,3], which has been observed 
ubiquitously in high-T c superconductors (HTSC's) for 
more than a decade. (See Ref. [3] for a review.) Sev- 
eral carefully-controlled experiments performed recently 
strongly supported the conclusion that the ZBCP's ob- 
served in them are due to such MS's [4]. 

Many other consequences of the MS's have been pre- 
dicted [5] , including a new contribution to Josephson tun- 
neling [6], a paramagnetic Meissner effect [7,8], and a 
low-temperature anomaly in the penetration depth [9], 
etc. The GMM is also a consequence of the MS's, which 
has not yet been looked for experimentally, perhaps be- 
cause no detailed analysis has been made on whether and 
how it can be observed. Thus here we perform such an 
analysis. 

An applied magnetic field B can have two simultane- 
ous effects on the surface MS's [1,7,9]: (i) Spin shift: 
The MS's are spin eigenstates. The field B can cause 
the MS's of one spin to shift above the Fermi surface, 
and those of the other spin to shift below. If the orbital 
shift (described below) can be neglected, then when the 
magnitude of the spin shifts grows past: (i) the width of 
the MS's peak in the DOS, (ii) the small non-zero ener- 
gies ~ Aq/Ef of the MS's, and (iii) the thermal energy 
fcgT, then a measurement of the total magnetic moment 
of the system should exhibit a saturation phenomenon of 
a magnitude proportional to the total number of MS's 
on the surface. This is referred to as GMM in [1]. (ii) 
Orbital shift: (including screcening current effects.) In 
the gauge in which the pair-potential order parameter A 
is real, this effect is due ot a vector potential A alone. 
The MS's acquire energy shifts proportional to their mo- 
menta k along A. At sufficiently low T all occupied MS's 
have the same sign of k, (in the WKBJ approximation 
and neglecting the spin shift,) implying a paramagnetic 
equilibrium current. Higashitani [7] thus proposed that 
this effect can account for the observed "paramagnetic 
Meissner effect" [8]. Another predicted consequence is 
a low temperature anomaly in the ab-plane penetration 
depth X a b which has been observed [9] . 

Both types of energy shifts are really present at the 
same time. Thus it is important to estimate their rel- 
ative magnitudes. We find that the conclusion depends 



crucially on the direction of B. Consider a thick single- 
crystal slab with a {110} surface, and with B applied 
parallel to the surface. If B is along [001], the screen- 
ing current is along [110], which is denoted as the y axis, 
with the x axis perpendicular to the surface at a; = 0. 
The spin shift has essentially the magnitude fioB. (/io is 
the Bohr magneton.) The orbital shift can be estimated 
using first order purturbation theory. The purtubation 
Hamiltonian is — (l/c)j- A(x). In the gauge in which A is 
real, A(x) = -BX ab exp(-x/\ ab )e y ~ A(0) = -B\ ab e y , 
because the MS's are localized within roughly one coher- 
ence length (£ a {, ~ 15A) from the surface, which is much 
smaller than X ab ~ 1500A (For the same reason the spin 
shift sees simply the applied field B.) Thus the orbital 
shift is ~ (e/m a bc)\k\B\ ab . The ratio between the or- 
bital and spin shifts is 2(m e /m b)(|fc|A a b). The mass ra- 
tio m e /m a b is probably less than unity by a factor larger 
than 0.2. [10] k ranges from — kp to kp, and kp should 
be somewhat less than w/2a, where a is the lattice pa- 
rameter in the ab plane. Thus this ratio is around 200, 
showing that with B along [001] the spin shift is negligi- 
ble in comparison with the orbital shift, and the previous 
analyses of the consequences of the orbital shift without 
considering the spin shift [7,9] is justified. 

Next let B be along [110]. The screening current is 
then along [001]. The same ratio in the simplest esti- 
mate is now 2(m e /m c )(k z \ c ). The London penetration- 
depth formula implies that m c oc A^. Thus this energy- 
shift ratio is reduced from the previous one by a factor 
(A q6 /A c ), which is - 1/50 for Hg-1201 [11]. More careful 
estimate, taking into account (i) the tight-binding nature 
in the c direction, (ii) k z ranges from —ir/c to ir/c with 

c/a ~ 3, reduces the energy-shift ratio to ^2. We shall 
see that this is quite sufficient to allow the GMM to be ob- 
served. On the other hand, by forming NS superlattices, 
it should be possible to increase A c by another factor of 
10 or larger, then one can even explore the regime where 
the orbital shift is negligible. Below, we first consider the 
effects of a spin shift alone, and then comment later on 
the effects of a, simultaneous orbital shift of a compara- 
ble magnitude. We shall discuss the spin magnetization 
first, and then the spin polarized tunneling conductance. 
We believe that the latter is a very promising way to see 
this effect. 

Consider first the spin magnetization M (which we de- 
ifne as the magnetic moment per unit area per Cu02 
plane). With 

M = - Mo [< ip\(x)^(x) >-< V{(a#|(a;) >], (1) 

where ips{x) is the field operator of spin-s electrons, 
its main contribution from the MS's follows easily from 
a perturbative treatment of the Bogoliubov-de Gennes 
(BdG) equations [1,3], which gives 

M(x)=nog(noB)f(x), (2) 
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E + iS lu- E + i5 
(3) 

k ^[\u k (x)\ 2 + \v k (x)\ 2 ], (4) 

where (3 = l/fc^T; u k (x) and ^fe(x) are the electron and 
hole components of the wave function of the MS of mo- 
mentum k along y; and S > takes into account a finite 
life time due to surface and/or bulk scatterings. In the 
limit S — ► + , g(E) — > tanh(/?_E/2). In general, 5 is a 
function of both P\xqB and /3(5, and is a monotonically 
decreasing function of T at any given B and <5. The 
largest value for g is unity, corresponding to [IqB >> {S 
and ksT}. Then the total magnetic moment per CuC>2 
plane, associated with the MS's on one {110} edge of 
the plane, obtained by integrating M(x) over x and y, 
is equal to (2L y / Xp)no- That is, for every Fermi wave- 
length on each {110} edge of a CuC>2 plane, there are 
two MS electrons contributing to the GMM. This is the 
maximum magnitude of the saturation phenomenon men- 
tioned earlier. To observe it directly, however, one needs 
to drastically increase the surface to volume ratio in the 
sample. The approach adopted in Ref. [9], where a sam- 
ple is irradiated with ions along [110], in order to create 
a large number of straight tracks, might offer some hope. 

Next we consider tunneling between a spin-polarized 
normal metal (N sp ), (i.e., a ferromagnetic metal or a half- 
metallic magnet [12],) and a DWSC (S d ) with a {110} 
surface. With B along [110] and the orbit shift neglected, 
the ZBCP should split into two peaks at non-zero volt- 
ages where the Fermi level of each spin species in N sp 
matches the shifted energy of the MS's in Sd of the same 
spin. The relative heights of these two peaks should de- 
pend directly on the spin polarization in N sp . Below we 
make these statements more quantitative by considering 
the effects of finite peak width and temperature. 

We assume, for simplicity, that both sides of the junc- 
tion have the same carrier type and density, and the same 
two-dimensional band dispersion (with HTSC's in mind). 
We also assume that the carriers are electrons with charge 
— e < and a gyromagnetic ratio 7 = —ge/2m e c with 
g = 2. Later we will comment on the effects of replacing 
these assumptions by more realistic ones, such as a three- 
dimensional band dispersion for N sp , and different carrier 
types and densities in the two sides of the junction. 
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The zero-field polarization in N sp is defined to be: 

P = (n T -n i )/(n T +n{) 



(5) 



where n s is the density of spin-s electrons at B = 0. 
We consider P > 0. The zero-field Fermi energies 
for spin-s electrons in N sp can then be expressed as 

1 + sP [s = 1(— 1) for spin 1(1)}, where 
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ergy Ep in Sd due to our simplifying assumptions. The 
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(Ep ^ + Ep^)/2, and is equal to the Fermi en- 



single-particle excitations in this system is governed by 
the BdG equations [13], 
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where H s (x) is equal to p 2 /2m + sfi^B — Ep\ for x < 
(i.e., in N sp , valid for Epl > [IqB only), and to 
p 2 /2m+sfi B(x)-Ep for x > (i.e., in Sd, B(x) ~ B for 
the MS's). The pair potential A vanishes in N sp , and is 
assumed to be x- independent in Sd- (Its self-consistency 
needs not be considered, for the MS's are topological.) e 
is the quasi-particle energy measured from the chemical 
potential, which is = Ep because we take the bottom of 
the conduction band in Sd to be zero. The bottom of 
the conduction band of each spin in N sp is then not zero, 

and has been absorbed in the definition of Ep \. In mo- 
mentum space, the d-wave pair potential A(kp) is taken 
to be A cos(2# — 2a), where 9 and a are the angles kp 
and the crystal a-axis make with the x-axis. We consider 
a = 7r/4. At a: = a 5- function barrier, H5(x), is as- 
sumed. All wave vectors are two-dimensional due to our 
assumptions. In N sp , the two-component quasi-particle 
wave function is, aside from a factor exp(iq y y), 



where V's.g, = ips exp(iq x x)/^\v qx | with ip^ 



(7) 



for 



a spin-up electron, and ip± 



for a spin-down 



hole, and v Qx — de/dq x is the qroup velocity. The 
first term in the right hand side of Eq. (7) expresses 
the incoming quasiparticle with spin-s and wave num- 
ber q x normal to the interface. The other two terms 
correspond to Andreev and ordinary reflections, respec- 
tively, with a s and b 8 their respective coefficients. Eq. 



(6) gives ti z q 2 x /2m = E^ s - s/j, B + se - ti z q 2 y /2m, and 

h 2 kl/2m = Ep- S + sfi B - se - h 2 q 2 /2m. The trans- 
mitted wave in Sd is a linear combination of outgoing 
waves (from the interface) that are solutions of the BdG 
equations for a bulk DWSC. Matching the wave function 
at the interface gives a s , b s , and other coefficients in the 
transmitted waves. 

The tunneling conductance G, (normalized to unity at 
Ep » eV » A when T, P, and B are all 0,) is cal- 
culated using the Blonder-Tinkham-Klapwijk formulism: 
[14,15] 

g =-£/o 00 ^/:5 2 # £ ^{E s =T,i 



f'(se + eV)[l + A s (e,<P) - B s (e,cj>)}}, 



(8) 



where 4> is the angle between q and the x axis, e(s, q, <f>) 
is given in the previous paragraph, qp = (2mEp) 1 / 2 /%, 
/(e) = l/[exp(e/fc B T) + 1], and A s = \a s \ 2 , B s = \b s \ 2 . 
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For numerical calculation, we take Ao/Ep = 0.08 as 
a typical value for HTSC's. The dimensionless barrier 

parameter is Z = j-\j2W^ ^ gi yes G as a 

function of V at Z = 5, for B = T = and for /j,oB = 
0.03A , ksT = (0, 0.2, and 0A)fi o B, for a non-magnetic 
N sp , P = 0. At B = 0, we get the ZBCP as observed 
in many experiments. At B ^ 0, the conductance peak 
splits into two peaks at eV = ±[ioB which correspond to 
the energy levels of the surface states of different spins in 
the SC. 

At [i B >> {ksT and the peaks' width), the spin- 
down surface states with energy e = —^qB are occupied 
with electrons, whereas the spin-up surface states at en- 
ergy +/j,oB are empty. As a voltage V > is applied 
between N sp and Sd, the chemical potential in N sp is 
lowered by eV. The spin-down electrons in the surface 
states in Sd will tunnel into N sp only when eV increases 
past noB, leading to a step-like increase in the tunneling 
current, or a conductance peak. If a negative voltage V 
is applied between N sp and S<j, the chemical potential 
in N sp is increased by e\V\. When it exceeds hqB, The 
spin-up electrons in N sp can then tunnel into the empty 
spin-up MS's in Sd, leading also to a peak in G. 

Figure 2 gives a similar plot for G at P = 0.5 (without 
the B = case, which does not change with P). The 
two conductance peaks at eV — ±/io-B now have different 
heights. The peak associated with tunneling of spin-down 
electrons is lower, because their Fermi velocity in N sp is 
smaller. (The DOS in two dimensions is a constant of 
energy, otherwise there would be another source for the 
peak-height difference.) 

Figure 3 gives G when N sp is fully polarized, P = 1, as 
found in some half metallic magnets [12] . In this case only 
one peak appears, which is associated with the tunneling 
of spin- up electrons. 

The absolute heights of these peaks have meaning, as 
G has been normalized. Note that the conductance peaks 
are higher for larger Z. At the same time, they become 
narrower. This is because, for higher barrier, the lifetimes 
of the surface states in Sd become longer, and there is a 
sharper resonance between the particles from N sp and the 
MS's in Sd- Therefore, to detect the GMM this way, it 
is preferable to work in the "tunneling limit", i.e., when 
the intcrfacial barrier is high. 

A typical value for A in HTSC's is about 16.5 meV. 
Then to reach the energy 0.03Ao, the magnetic field 
needs to be around 8.6 Tesla. In addition, 0.03Ao in 
temperature is about 5.7K. To measure such an energy 
shift, the experiment should be performed at liquid he- 
lium or lower temperatures. 

If an orbital shift of the same order as the spin shift is 
present in the system, then each peak will become wider 
and approximately rectangularly shaped. But the effects 
predicted here are still observable, even if the maximum 
orbital shift is, say, three times larger than the spin shift. 
Correction to the WKBJ approximation has been ne- 
glected, otherwise [i B might have to be larger to see 



this effect. (Accurate estimate of this correction is diffi- 
cult, but it is of the order Aq/Ef.) 

If the carrier density in N sp is different from that in 
Sd, or if the band dispersion relations of the two sides 
are different, then the main effects, as far as we can see, 
are: (i) a change in the absolute magnitude of the tun- 
neling conductance which does not affect the normalized 
conductance, (ii) a stronger P dependence due to the en- 
ergy dependence of the DOS in N sp (without changing 
the limiting behavior at P = 1), and (iii) a change in the 
ordinary reflection coefficient at the interface, which can 
be simulated with an effective Z. Even when the carriers 
in Sd are holes, and those in N sp are electrons, we still 
find no essential change in our predictions. 

In summary, we have analyzed in detail the giant mag- 
netic moment that can result from the midgap states 
on the {110} surface of a rf-wave superconductor. With 
high-T c superconductors most-likely having d-wave pair- 
ing, this predicted giant magnetic moment should be ob- 
servable either directly in samples with a high concentra- 
tion of internal {110} surfaces generated by ion irradia- 
tion, or, with better promise, by measuring the tunneling 
conductance between a spin-polarized normal metal and 
a high-T c superconductor with a {110} surface. At a 
large enough external magnetic field applied along [110], 
and low enough temperature, the ZBCP is shown to split 
into two peaks, with their relative heights determined by 
the polarization of the normal metal, because each peak 
is associated with a single spin. Observing these pre- 
dictions can confirm the existence of the midgap states, 
unconventional pairing, and remove any remaining doubt 
that the observed zero-bias conductance peak might not 
be due to the midgap states. 
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Figure Captions 

Fig. 1 Normalized tunneling conductance G between an 
unpolarized (P — 0) normal metal and a d-wave 
superconductor with a {110} surface, as a function 
of voltage V, for B = 0, T = 0, and B = 0.03A /mo 
at three values of T. The interfacial barrier parame- 
ter Z = 5. Only the contributions from the midgap 
states are included. 

Fig. 2 Similar plot for G except that the normal metal is 
50% spin-polarized (P = 0.5). Only the B ^ 
case is plotted since the B = case is practically 
independent of P. 

Fig. 3 Same as Fig. 2 except that the normal metal is 
100% spin-polarized (P = 1). 
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